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Abstract. We review and study a one-parameter family of functional transformations, 
denoted by (5''''')/3gE, which, in the case /3 < 0, provides a path realization of bridges 
associated to the family of diffusion processes enjoying the time inversion property. This 
family includes the Brownian motion, Bessel processes with a positive dimension and their 
conservative ft-transforms. By means of these transformations, we derive an explicit and 
simple expression which relates the law of the boundary crossing times for these diffusions 
over a given function / to those over the image of / by the mapping S'"'', for some fixed 
/3 £ R. We give some new examples of boundary crossing problems for the Brownian 
motion and the family of Bessel processes. We also provide, in the Brownian case, an 
interpretation of the results obtained by the standard method of images and establish 
connections between the exact asymptotics for large time of the densities corresponding 
to various curves of each family. 



Let X := {Xt,t > 0) be a E'-valued, where E = M or M_|_ = [0, oo), conservative, in 
the sense that it has an infinite hfe-time, 2-self-similar homogenous diffusion enjoying the 
time inversion property in the sense of Shiga and Watanabe [39]. We denote by {T'x)xeE 
the family of probabihty measures of X which act on C(M+,i?), the space of continuous 
functions from into E, such that Pa;(^o = x) = 1. The standard notation J^t is used 
for the (T-algebra generated by the process X up to time t and we write simply = T^o- 
The above assumptions mean that, for any x G E and c > 0, 

the law of {c~^X(,2f,t > 0;Pca;) is P^^ (2-self-similarity) 

and the process 



is a homogeneous conservative diffusion. It is well known that the class of processes 
considered therein consists of Brownian motion and Bessel processes of positive dimensions, 
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see Watanabe [H]. Next, let / G C{R+,E), such that /(O) / Xq and set 

= mf{s>0; = f{s)} 

with the usual convention that inf{0} = +00. The study of the distribution of the stopping 
time is known as the boundary crossing or first passage problem. Unfortunately, the 
explicit determination of these distributions is only attainable for a few specific functions. 
In this paper, we suggest a new method which allows to relate, in a simple and explicit 
manner, the law of T-f with the ones of the family of stopping times (T-^*'^')^gR with 
:= 5{/3)(/) and 

(2) S(^) : C{R+,E) ^ C([0,C(^)),S) 

where 



1 + 




if /? < 0, 
otherwise. 

The results extend to /i-transforms of the prescribed processes leading to conservative 
diffusions. 

In the case /? < 0, our methodology has a simple description. Indeed, as a consequence 
of the time inversion property. Pitman and Yor [34J showed that the law of the process 
{S^f^\X)t,0 < t < C^^^) corresponds to the law of a bridge associated to X. Thus, one 
may relate, by means of a deterministic time change, the law of with the first crossing 
time of the curve /(^) by the bridge {S^^\X)t,0 < t < C^^^). Finally, our identities are 
then obtained by using the construction of the law of bridges of Markov processes as a 
Doob /i-transform of their laws, see Fitzsimmons et al. (16J. We also show that similar 
devices can be implemented in the case (3 > 0. 

The motivations for investigating the boundary crossing problem are of both practical 
and theoretical importance. Indeed, on the one hand, such studies were originally mo- 
tivated by their connections to sequential analysis, see e.g. Robbins and Siegmund [37] 
and Smirnov-Kolmogorov test, see e.g. Lerche [26] • On the other hand, this problem has 
found many applications in several fields of sciences, such as mathematical physics [27], 
neurology [23], epidemiology [28] and mathematical finance [5] and |30j . 

Amongst general results, in the Brownian setting, Strassen [IQ] proved that if / is 

continuously differentiable then the distribution of T-^ is absolutely continuous with a 

continuous probability density function. Moreover, for elementary curves, the density is 

known explicitly for linear, square root and parabolic functions. For these curves, specific 

techniques proved to be efficient and we refer to Table [1] in Section 2 for a description of 

these cases. Besides, the distribution is known for some concave curves solving implicit 

equations obtained by the celebrated standard method of images, see Lerche [26]. For 
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some recent investigations, we refer to Potzelberger and Wang [35j and Borovkov and 
Novikov [7] for numerical approximations of the density, Peskir [32], [33j for the study of 
the small time behavior of the density and Kendall et al. [21] for statistical applications. 

The remaining part of the paper is organized as follows. In the next Section we recall 
some recent results regarding diffusions which enjoy the time inversion property. Section 
[3] is devoted to the statement of our main results and their proofs. Section H] is concerned 
with a detailed study of the Brownian motion case. In particular, we present several 
new explicit examples of the boundary crossing problem. We also show how our results 
translate and agree with the standard method of images. Finally, in Section [H we treat the 
example of Bessel processes and characterize the distribution of hitting times of straight 
lines. We also mention that, in the Brownian motion case, the boundary crossing identity 
([7]), stated below, has been published without proofs in the note [I]. 



2. Preliminaries 

Let us recall that E is either R+ or R and X := (Xf, t > 0) is a 2-self-similar conservative 
homogenous diffusion enjoying the time inversion property. In recent papers, Gallardo 
and Yor [18] and Lawi [24] characterized the class of self-similar Markov processes (which 
might have cadlag paths) satisfying the time inversion property in terms of their semi- 
groups when the latter are assumed to be absolutely continuous and twice differentiable. 
More precisely, they showed that if we write Fx{Xt G dy) = pt{x,y)dy the transitions 
densities have the following form 

pt(x,y) = ^$(^) (^^) e-'V, t>0,x,yeE, 

where the reals and c are related to the function ^ : E ^ as follows. 

(1) If = M, then X is a Brownian motion, <I?(y) = and necessarily u = —1/2 and 
c= l/\/2^. 

(2) If E = M"*" then X is Bessel process of dimension 6 for some 

and ^{y) = y^^Ip{y) where ly is the modified Bessel function of the first kind of 
index v which admits the following power series representation, see e.g. [221 Section 
5.2], 



n=0 ^ ' 

Moreover, it is well known that, in this case, the state is 

- an entrance boundary if (5 > 2, 

- a reflecting boundary if < 5 < 2. 
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Next, recalling that, for every fixed y £ E, the mapping x ^ ^(yx) is ^-excessive, see 
e.g. [6], we can define a new family of probability measures as follows 

(3) dpr,^.=^.-i»"dp.,^., *>o^ 

Note, from the definition of that = P^,. Moreover, for any y £ E, {^^^)x(^e is 
precisely the family of probability measures of the process i{X) when X starts at y. We 
also recall that i{X) is conservative and satisfies the time inversion property ([T]), see |18) . 
It follows that i{X) has an absolutely continuous semi-group with densities given, for any 
t > and x,a £ E, by 

(4) p[y\x,a) = ^i(4$r^U^V'^^\-M*^^+(-^+<^^)A). 



^/i^{yx) V t J 
Note, in particular, the following: 

(1) If = R then i{X) is a Brownian motion with drift y. 

(2) li E = M"*" then i{X) is the so-called Bessel process in the wide sense introduced 
by Watanabe [4lj. 

Now, from the definition of the family of mappings {S^^'^)f^>Q, see ([2]), we observe that 
for any / £ C(M+,M) and a,/3 G M, we have 5'W(/) = / and S^°''> o 5'(^) = 5^"+^) and 
so {S^^^)j3>Q is a semi-group on C(M^,M). Moreover, the family of mappings (5'^^^)/3<o 
naturally appears in the process of construction of bridges associated to X. To explain the 
connection, for any fixed T > and x £ E, let us denote by (P^^^ ^, z £ E) a regular version 

of the family of conditional probability distributions (f^x\.\Xt = z), z £ Ej . That is, for 



(v) 

some z £ E, P^^ j, is the law of the bridge of length T associated to X, between x and z. 

This, according to Fitzsimmons and al. [16j, can be obtained from pi?'^ as the following 
Doob /i-transform 



(5) Fil^ = ^T^IPi^^ on^.,.<T. 

p)f [X, z) 



Fitzsimmons [17J observed that these laws remain invariant under changes of probability 
measures of type ([3]). Finally, we recall that Pitman and Yor [341 Theorem 5.8], see also 
[181 Theorem 1], showed, by means of the time inversion property that, for any x, z £ E, 
the processes 

(6) {Xu,0< u < T;F^^}r^^^} and {S^-^/'^\x)u,0 < u < T;Fi'^} 

have the same law. Note that, in the Brownian case, ([6]) can also be seen from the unique 
decomposition, as a semi-martingale, in its own filtration, of the Brownian bridge of length 



4 > 0, between and 0. Indeed, for any < t < T, we can write 



/3 



Jo T-s Jo T~s 

where X is a Brownian motion on [0, C''"'^^] with respect to the filtration generated by 
(5("^)(X)t,0 <t<T). Thus, we have 

S^^y^){X), = X,- f ^^-^^ds, t < T. 
Jo T-s 

which coincides with the canonical decomposition of the standard Brownian bridge. 

3. Main results and proofs 

We keep the notation and setting of the previous section and assume, throughout this 
section, that /? is some fixed real. We proceed by pointing out that the mapping can 
be defined similarly on the space of probability measures. For instance, in the absolutely 
continuous case, we associate to fi{dt) = h{t)dt the image S^^\fi)idt) = h^^\t)dt where 
we recall that h^^\t) := S^^\h){t). We are now ready to state our main result. 

Theorem 3.1. Let f £ C(M+,£'). Then, for any x, y G E such that /(O) / x, and 
t < C^^\ we have 

X s^'^'^ (F^y\Tf e dt) 

The particular case y = yields 

(7) (Tf''' G dt) = (1 + /3f)-2e-§^S^+^' X 5^^) (p, (t^ G dt 



Theorem 13.11 simplifies when the focus is on straight lines. Indeed, if we consider a 
constant function f = a where a 7^ then, with (3 = b/a, for some 6 G M, we have 
= a + bt, t < C^*/"). Note that if 6 < then the support of r"+^- is (0, -b/a). 
The previous result reads as follows. 

Corollary 3.2. For any fixed a,b,x £ E such that a ^ x and all t < cC*/"), we have 

P(^)(r--Gdt) = (l+6t/ar^iM^exp-^fa + 6t+^-:^) 
V / ^{ya) 2V a + bt a J 

X (^p^s')(r" G dt)) . 

In particular, when y = 0, we obtain 

(r'^+^- G dt) = (l + bt/ay-^exp-^ia + bt-x^a) 
X (P^(r" G dt)) . 
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Before proving Theorem 13.11 we need to prepare two intermediate results. To this end, 
let us denote by H^^^'^^ the first time when S^^\X) crosses / i.e. 

H^-M = inf {o < 5 < C^^); S^^\x)s = /(«)} . 

The aim of the next result is to relate the stopping times and T^^'^^ . 

Lemma 3.3. The identities 



1 + /?r/<« 1 - f5Hi-M 

hold almost surely. In particular, we have ^H^~^'-^^ < = It-^'*^' < C^'^^j- 

Proof. From the definition of S^^\X) and by using a deterministic time change, we get 

H^-M = m{o<,<c'-«;X^ = ^} 



1 + PT'"> 

The second identity is obtained in a similar way. The last statement follows then by 



observing that, for any /3 G M, we have C*-^^ = ^ "^grf-zS) ™ limiting sense. □ 



The image by the mapping 5^'^^ of any homogeneous Markov process is clearly a non- 
homogeneous Markov process. However, the time inversion property allows to connect the 
law of X and that of S'(^)(X) via a simple time-space /i-transform. To be more precise, 
we state the following result. 

Lemma 3.4. Let x,y e E and (3 e R. Under Fi^\ the process X^l^^ = S^^\X), defined 
on [0,C^^)), is a non-homogeneous strong Markov process. Its law, denoted by Qx'^ , is 
absolutely continuous with respect to P^f^ with Radon-Nikodym derivative Mx'^ {t, Xt), t < 
C^^\ given by 



g 2 l+0t~'~ 2 l+/3t 2 



In particular, for y = 0, we have 



M^''^{t,Xt) = (1 + /3t)"''"^e2^"2^'. 



Proof. We start with the case /? < 0. Using successively the identities ([6]) and ([5]), we 
obtain, for any bounded measurable functional F and t < C^^\ that 



E^y^[F{S^''Hx)s,s<t)] = E^y^[F{X„s<t)\X_^ 



in 



E(^) \My^^{t,Xt)F{Xs,s<t) 



F{Xs,s< t) 



where the last line follows from (jj]). Using the fact that, for any c,d £ E, 



1 as 



y ^ 0, we get the stated result for y = 0. Next, we treat the case (3 > and assume, first, 
that X is a Bessel process of dimension 5 > and y = 0. Using the Ito formula, we see 
that the process {Mx'^ {t, Xt),t > 0) is a P2.-local martingale. We now show that it is a 
true martingale. Recall that the squared Bessel process is the unique solution to the 
stochastic differential equation (for short sde) 



(8) 



dX^ = 2XtdBt + 6dt, X?. 



where i? is a Brownian motion, see |36j . Furthermore, we need to introduce the processes 
Y and Y which are defined, for a fixed t > 0, by 

Yt = {1 + pt)Yt = S^fHx)t. 

It follows from ([8]) and by performing a deterministic time change, that Y solves the sde 



dY/ 



dt + 



2Yt 



(l + /3t)2 ' 1 + pt 



-djt 



where (7f,t > 0) is a Brownian motion with respect to the filtration generated by the 
process Y. Similarly, Y'^ satisfies 



(9) 



((i + /3t)2y/ 



d 



2Yt djt + 



l + Pt' 



Ytdt + 5dt. 



Now, from ([8]) and by Girsanov theorem, the law of the unique solution to Q is ob- 
tained from the law of X"^ by a change of probability measure using the local martingale 
{Mx'^ {t, Xt),t > 0). The proof of the claim, in the case y = 0, is completed by invoking the 
conservativeness property of y^, which implies that {Mx'^ {t, Xt),t > 0) is a martingale. 



Finally, to recover the case y 7^ 0, we use ^ which gives 



= ^w(i+w) ,-ijf;,My«.x.)dF.i.. 

This completes the proof in the Bessel case. The same arguments work for the Brownian 
case and this completes the proof. □ 

We shall now return back to the proof of Theorem l3.11 Lemmas l3.3l and l3.41 when combined 
with the optional stopping theorem, allow us to write, for any A > 0, 



(/3) 



I 



{T/"''<C('3)} 



1_/3h(-/3,/) J 



{jy(-/3./)<C(-«} 



-X — 

e 1- 



-Am 



(l + /3n)3 



where we have performed the change of variables t/{l — (3t) = u. Simplifying the above 
expression and by the injectivity of the Laplace transform, we get our first assertion. The 
second statement follows by letting y — > 0. 



4. Brownian motion 

In this paragraph, we take E = M. Thus, X is a Brownian motion. Thanks to the 
homogeneity property, in this case, it is clearly enough to consider the case Xq = and we 
simply write P for Pq. To simplify the discussion, except in Subsection l4.3l we assume that 
/ G C"'^(]R4-, E), with /(O) / 0, which implies that the studied distributions are absolutely 
continuous with respect to the Lebesgue measure with continuous densities, see [40]. We 
write then P (T-^ G dt) = p-^ {t) dt and read from Theorem 13.11 the identity 

^'"«-(TT^^-™""'"V(t^).*<c-. 



Observe, when f3 > 0, the asymptotic behavior 

(11) pf''\t) ~ (/3i)-3/2e'^i^^''''W/(l//3) 



as t 



oo, 



where h{t) ~ g{t) as t — > means that h{t)/g{t) ^ 1 as t — > for some C, € [0,oo]. It 
seems natural to examine the images of some curves by the mapping 

We gathered m 

Table [U the images of the most studied curves. The fact that S^'^'i preserves straight lines 



/ 
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Table 1 . Image by S^^^ of / and the corresponding references where the 
distribution of is studied with a, 6 G M, 6i > —6^/4 and b = bie~°-^^ . 



is well known, see for instance [34J. More generally, we observe that if gp^a{-) = (1 + P-)^) 
for some reals a and /3, then we have ff^^^ = g-p^i-a- In particular, if we take f = a 
and /? = 6/a, for some reals a and 6, then f^^^'^\t) = a + bt, t < (,^^1°^ . An immediate 
application of Corollarv 13.21 combined with 

^{T" G dt) = Jj— e-'idt, t > 0, 

yields then the following well known Bachelier-Levy formula 



4.1. Some new examples. We now consider the boundary crossing problem associated 
to two families of curves consisting of the square root of second order polynomials and the 
reciprocal of affine functions. 

4.1.1. First, we consider the distribution of the stopping time 

T(^i,X2) = inf |s > 0; X, = aV(l + Ais) (1 + Ass)} 

where a and Ai < A2 are some fixed reals. We do not consider the case Ai = A2 which can 
be studied in a elementary way, with the extra cost of making use of the strong Markov 
property when Ai < 0. First, consider the case A2 = and, to simplify the notation, set 
Ai = A and Ta^'^^ = Ta^\ This is the setting of the classical example studied by Breiman 



in [S]. T^^^ is related to the hitting time of a constant level by an Ornstein-Uhlenbeck 
process and we refer to Alili et al. [2] for a recent survey on this topic. That is, with 



and Ha = inf{s > 0; Us = a}, we have 
(12) rW^=^A-i(e^^'"-l^ 



where == stands for the identity in distribution. By symmetry, it is enough to consider 
the case where a is positive. We proceed by recalling that the distribution of Ha is given, 
see for instance [2], by 

1 D (0) 

P ^Ha e dt) = -iAe-^»V4 g '^"--"^^ ^ e-^^".-^^*/^ t > 0, 

where we used the notation L'^/^\(6) = ^%7^|i/=i/„b and {i^j,b)j>o stands for the ordered 
sequence of the positive zeros of the parabolic function v Dy{h). By means of the 
identity (fT2]) . we get that 

p (rW G dt) = Y^P(^^a e rf-) Li i > 0. 

Next, we assume that Ai < A2. Thus, the support of T^^^'^'^'' is [0, C^"^^^) if Ai is positive 
and is the positive real line otherwise. We have 

(Vl + (A2-Ai)-) = x/(l + A2-)(l + Ar)- 
We are now ready to use Theorem 13.11 and write 

4.1.2. We are now interested in computing the distribution of the stopping time T^^'^'' 
defined by 



r-^i„,{o<.<c<«;X,^^} 



where [3 is some real. To this end, we recall that Groeneboom [19j has computed the 
density of with h{t) = 1 + I3H'^ as follows 

P(r^ e dt) = 2{p\fe-l''''' f; Mfk + 2c(3^^_,^, ^ ^ 

where {zk)k>o is the decreasing sequence of negative zeros of the Airy function, see e.g. [25], 

and we have set c = (2/3^)~3. Next, by means of the Cameron-Martin formula, we obtain, 
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with h{t) = (1 - Ptf = h{t) - 2f3t, 

FiT'^ G dt) = 2(/32c)2e2,3-2/32<i+§,3^t^-/3t) y M^k + 2c0^^_^^, ^ ^ ^ 

Finally, since h^^^ = S^^\h), we obtain, by applying Theorem 13 -H that 
We complete the example by stating the following asymptotic result 



P(r^"" G dt) ~ f 2(/32c)2e2'^ f; M^_t^e-T] dt as t 



Ai'(zk) 



oo 



which holds provided that /3 > 0. 

4.2. Interpretation of the mapping via the method of images. We aim to 

describe the impact of our methodology to the so-called standard method of images. To 
this end, let us assume that Xq = and set 

h{t, x)dx = f{t^ > t, Xt G dx) . 

We need to impose some conditions on the the curve / in order to be able to apply the 
standard method of images. That is, we assume that 

- / is infinitely often continuously differentiable, 

- f{t)/t is monotone decreasing, 

- / is concave. 

Note that these properties are also satisfied by /^^^ when P > 0. The function h is 
characterized as being the unique solution to the heat equation, see Lerche [26l Chap. I.l], 
dh 1 d'^h 

m=2d^ ^) e X M; X < fit)}, 

with the boundary conditions 

h{tj{t)) = 0, h{0,.)=6o{.) on ]-oo,/(0+)] 

where 5q stands for the Dirac function at and /(O^) = limj^o /(O- 'T^^ standard method 
of images assumes that h is known, whilst / is unknown, and is given by 

1/2 /^OO j _ .2 

(13) h{t,x) = ^={e-^ - e~-^F{ds) 

where F{ds) is some positive c-finite measure satisfying e~~ F{ds) < 00, for all e > 0. 
In [26], it is shown that if / is the unique root of the equation h{t, x) = in the unknown 
X for a fixed t >0, then we have 

Tfedt]=^= se-^'^^Fids), t>0. 
' V27rt3 Jo 
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Note that the imphcit equation h{t, x) = may be written as 



(15) / e-^+'fF((is) = 1. 

Jo 

With F{ds) replaced by F{ds)e~^^ , fi > 0, the unique solution to is f{t) + fit for 
a fixed positive real t and this is easily checked by the Cameron-Martin formula. In the 
same spirit, the identity pO|) has the following interpretation. 

Proposition 4.1. > then we have the following representation 

/ Jo 

Furthermore, the a -finite measure corresponding to the curve 

f(P) isF^^\ds) := e-P''/^F{ds) 

and 



Ti'" >t,Xt€dx]= [e-2t - I e~^-2t^e-'' 2F{ds] 



Proof. Noting, from that 



T/'^edt) = I se-'^^'^FiPXd. 



dt («-/W(t))2 



V27rt3 Jo 

and by using ([TO]) we obtain that 



IS 



se 



'-F(^\ds) = _±_-e-^T^^'''WV f- 



1 /""^ (.-/('^'W)^ 



1 — e ai-f/st/'^'W^ / se 2t/{i+/3t) " 



se 2t e 2 F{ds). 

Jo 



\/27rfi Jo 

This proves the first assertion and suggests, but doest not prove, that the cr-finite measure 
corresponding to the curve 

fW) is F^^\ds). Next, with 

h(0)^t,x) = (e-i - r e--^e-P''^I^F{ds) 

v27rt V Jo 

we aim to solve h^^\t^x) = for each fixed t > 0. Now, setting x = f{t) and replacing x 

and i, respectively, by f{t) and t/{l+(3t) in 1^, we find that e"^+^^^^e-^^'/2^((is) = 
1. It follows that f^^\t) solves the implicit equation h^^\t,x) = 0, t > 0, which finishes 
the proof. □ 
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4.3. Large asymptotics for the density of T^"''. Following Anderson and Pitt [3], we 
consider the asymptotic of the density of the distribution of T-^'*^' as time gets close to 
(^^^\ To start with, we exceptionally assume in this subsection that / € C^((0, cxd),M+). 
So the distribution of T-^ has a continuous density with respect to the Lebesgue measure 
on (0,00). Assume that the distribution of T-^ is defective. That is r = ¥{T'^ < 00) < 1 
and, in this case, we say that / is transient. By the classical Kolmogorov-Erdos-Petrovski 
theorem, see [H], we know that if t-^/'^f{t) is increasing for sufficiently large t, then / is 
transient if and only if the integral test 



(16) / r^'^f{t)e-f"^'^'^'dt< 



00 

holds. Clearly if /? > and < f3f{l/(3) < 00 then /(^) does not satisfy ([ISD but the 
general formula (jlip holds. Now, we consider the more interesting case /? < and examine 
the asymptotic of the density of the distribution of T-^"'' as t — > — We need to work 
under the following three conditions borrowed from [3]: 

- / is increasing, concave, twice differentiable on (0, 00) and of regular variations at 
00 with index a E [1/2, 1), 

- f{t)/Vi is monotonic increasing at 00, and f{t)/t is convex and decreases to for 
sufficiently large t, 

- There exist positive constants c < 1 and c' such that we have the inequalities 
tf'{t) < cf{t) and \t'^f"{t)\ < c'f{t) for a sufficiently large enough t. 

The behavior at 00 imposed in the above conditions is granted in the examples where / 
behaves like 

fit) = Cr(logt)''(loglogt)^(logloglogt)'^ 

with 1/2 < a < 1, for large t. It is clear that these conditions are not always preserved by 
the family S^^\ Equipped with this, we are now ready to state the following result. 

Proposition 4.2. Assume that f is transient and satisfies the above conditions. Then, 
for any /? < 0, we have 

(t) ~ ^ ^ (1 - r) f{(3, t) ast^ C^^^ 

where 

mt) = f(T-^]--r^J'^ ' 



1 + I3tj l + (3f \1 + I3t 
Proof. We read from Theorem 1 that 

(',^\ ft.\ /I N /(*) ~ -f^(t)/2t , 

(17) f(^) ^ [l-r) /^^3/2 e ast^oo. 



Moreover, as t ^ 00, fit) /t I and hence 
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The proof is then completed by combining (|17p with (jlOp . 



□ 



Remark 4.3. In the defective case, the distribution of the last crossing time = sup{s > 
0;Xs = f{s)} is shown to have an atom at and its asymptotic as t ^ is determined 
in [lO]. Similar questions can be treated for using this 

method. 



5. Bessel processes and straight lines 



We investigate here the case when X is a Bessel process of dimension 5 > and we refer 
to Revuz and Yor |36^ Chap. XI] for a concise treatment of these processes. For 6 > 2, 
and X > 0, the Bessel process of dimension 6 is the unique solution to 

6-1 ds 



Xt = Bt + x + 







t > 0, 



where i? is a standard Brownian motion. For < 6 < 2, X is defined as the square root 
of the unique non- negative solution of ([8]). We recall that is an entrance boundary when 
5 > 2 and a reflecting boundary otherwise. We denote by PJJ (resp. E^) the law (resp. the 
expectation operator) of a Bessel process of index u = 6/2 — 1, starting at x > 0. The 
semi-group of X'^ is characterized by 



(18) 



-xxf 



(1 + 2\t) 



-S/2 



e i+2^t . 



A,t > 0. 



The densities of the semi-group of X, with respect to the Lebesgue measure, are 
Pt{x,y) = 7 - e 2t 1 t,x,y>0. 



For a given / G C(R+,E) let us keep the notation used in the introduction and Section [3l 
Observe, that if y = then Theorem 13.11 reads 



G dt 



(1 +/3t)'^-2g-§i^^ + fx2^(;3) (z dt)'^ 



for all t < C^l^\ We end up our discussion by computing the distribution of the hitting 
time by X of a straight line a + b- where a > and 6 is a real. We keep the notation used in 
Corollary 13.21 and the reader is reminded about observations preceding its statement. To 
the best of our knowledge, the problem of the determination of the distribution of T°'+^-^ 
which was raised in [34j for the case a ^ 0, remained open. Recall that the law of T", 
which corresponds to 6 = 0, is characterized by 



(19) 



_ A rpa 

e 2 ^ 



'K^{aX) ■ 



X < a, 



X > a, 
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for any A > 0, where Ky is the modified Bessel functions of the second kind of index see 
for instance Borodin and Salminen [6]. Observe that when x > a the distribution of T*^ is 
defective and IPJ5(r" < oo) = [ajx)^^ . It is also weh known that if x < a then we have 



(20) 



CxD ~y • J ( ' X\ 



where {ju,k)k>i is the ordered sequence of positive zeros of the Bessel function of the first 
kind Jjy(.), see [6]. Furthermore, if a = and 6 > then T^' and l/sup{s > 0;Xs = b} 
have the distribution which was determined in [34] by making use of the time inversion 
property. Now, we are ready to state the following result. 

Theorem 5.1. For < x < a and 5 G R, we have for any t < (^(^/°) 



e2a 



(1 + 1*)'^+' t^i a2-J,+i(j,,,) 
For any x > and a,b > 0, we have 



dt. 



e 2 



e2a 



a- Jo I4V2au)Pb/2aK^b,U)au, X 



< a, 



^ a-" h KA,^au)Phl2a\^b,U) au, X ^ a. 



where = {\^ + h'^) /2, A G M. 

Proof. The first statement results from a combination of Corollarv 13.21 and relation (j20p . 
Next, using Lemmae 13.31 and 13.41 with j3 = —b/a, we can write 



e 2 



rpa—b- ^ 



,{b/a,a) 



2 T+It" / ^ _|_ :i J^a 



a 



-S/2 b_ 

2a \ -|^_|_ o_ rpa 



\2t" 



A(a2_x2) 



6 2 







a 



du 



where the last line follows from (1181). It remains to use (1191) to conclude. 



□ 



Remark 5.2. The process (1^ := + bt,t > 0) is a non-homogeneous Markov process 
and solves the sde 

r. 6-1 ds 
Yt = Bt + I + bt, t> 0. 



Ys-bs 
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This is to be distinguished from a Bessel process with a "naive" drift b introduced in 
and defined as a solution to 

Zt = Bt + ^ f ^ + bt, t> 0. 

Remark 5.3. The Bessel process of dimension 5 = 1 is, in fact, the reflected Brownian 
motion. It follows that the associated first hitting times can be interpreted as double 
barrier hitting times. That is, with / as above, the time when a Brownian motion B hits 
one of the curves x ± /(•), i.e. inf{s > 0; Bg = x ± f{s)}. 

6. Concluding remarks and some comments 

6.1. It is plain that the results of Theorem l3.1l can be readily extended to any /i-transform 
of the process X, but we need to take care of the life-times of the involved processes in 
the Bessel case. In the Brownian setting, one gets similar results for the process = 
Xt + yet, t > 0, where e is an independent symmetric Bernoulli random variable taking 
values in {—1,1}. Observe that i(X^) is a strong Markov process. However, because the 
latter starts at the random point ye, X^^") does not satisfy the time inversion property ([T]). 

6.2. Assuming that X is a 2-self-similar strong Markov process with cadlag paths (i.e. with 
possible jumps) enjoying the time inversion property, then Theorem 5.8 in [34J or Theorem 
1 in [TH] allow to extend Lemma 13.31 and Lemma 13.41 and an analogue of Theorem 13.11 can 
be stated. However, we did not succeed to construct examples of M+-valued processes en- 
joying the time inversion property other than Bessel processes in the usual or wide sense. 
This is the reason why the setting is restricted to the continuous one. We refer to |31j 
where the second author characterizes, through its Mellin transform, the law of the first 
passage time above the square root boundary for spectrally negative positive 2-self-similar 
Markov processes. 

6.3. In [12j, Durbin considered the studied boundary crossing problem for continuous 
gaussian processes and showed that, in the absolute-continuous case, the problem reduces 
to the computation of a conditional expectation. In particular, for the Brownian motion, 
if / is continuously differentiable and /(O) 7^ then 

1 £m. 



P (t^ G dt 



2t 



/27rt 



'h{t)dt 



where 



h{t) = hm-^E \Bs - f{s);Tf > s\Bt = f{t) 

s/tt — S I ' 



t > 0. 



There seems to be no known way to compute the function h. It is shown in [13] that this 
method is in agreement with the standard method of images. With the obvious notation, 
it is clearly possible to relate h^^^ to h. 
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6.4. We learned from Kendall [20] an intuitive interpretation of Durbin's expression in- 
volving the family of local times of X at / denoted by (L^ ~ ,t > 0). Indeed, observe 
that Durbin's formula can be rewritten as F {T^^ € dt) = /i(t)E[dLf "^], t > 0. The 
latter when integrated over [0,a], yields F {T^ < a) = E \f^h{s)dL^-n. Such a de- 
composition is not unique and many can be constructed from the above one. A natural 
non trivial one to consider is motivated by the following observation found in [20]. If 
(/ : M+ X M+ ^ M+ solves E g{s,a)dLf"^ \Xt = f{t) = 1, for < t < a, then clearly 

E g{s, a)dL^~^ = F (T^ < a) . However, it is not clear how to express g in terms of 
h. W. Kendall told the first author that the above interpretation could be interesting for 
further investigations for the studied boundary crossing problem. 

6.5. In the Brownian case, it is interesting to analyze the impact of our identities on 
some integral equations satisfied by the studied densities. However, some of them lead to 
obvious facts explainable by change of variables. As an example, we observe that if / is 
positive and does not vanish then X, when started at a; > /(O), must hit / before reaching 
0. The strong Markov property then gives 

^e-i= fFATUdr)^M 
V2^ Jo A/27r(t 

This is easily shown to be in accordance with the result stated in Theorem 13.11 For the 
above and other classical integral equations we refer to [12], [15], [26] and also to [32j for 
some more recent ones. 

6.6. The technics used in the example treated in Subsection 14.11 Brownian motion, 
can be applied to Bessel processes and square root curves. The required results, for that 
end, can be found in Delong [llj or in Yor |42j . This law is connected via a deterministic 
time change to the one of the first passage time to a fixed level by the radial norm of 
a (5-dimensional Ornstein-Uhlenbeck process, where 5 is some positive integer. Due to 
the stationarity property, the law of the first passage times can be expressed as infinite 
convolutions of mixture of exponential distributions, see Kent [22]. 
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